arXiv:l504.06228v2 [math-ph] 25 Nov 2015 


Symmetry, Integrability and Geometry: Methods and Applications 


SIGMA 11 (2015), 096, 23 pages 


Harmonic Oscillator on the SO (2, 2) Hyperboloid* 

Davit R. PETROSYAN t and George S. POGOSYAN ^ 

Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 

Dubna, Moscow Region, 141980, Russia 
E-mail: petrosyan@theor.jinr.ru 

^ Departamento de Matematicas, GUGEI, Universidad de Guadalajara, 

Guadalajara, Jalisco, Mexico 
E-mail: george.pogosyan@cucei.udg.mx 

§ International Genter for Advanced Studies, Yerevan State University, 

A. Manoogian 1, Yerevan, 0025, Armenia 
E-mail: pogosyan@ysu.am 

Received April 24, 2015, in final form November 20, 2015; Published online November 25, 2015 
http://dx.doi.org/10.3842/SIGMA.2015.096 

Abstract. In the present work the classical problem of harmonic oscillator in the hyperbolic 
space iJ|: z^+zl — z^ — zj, = has been completely solved in framework of Hamilton-Jacobi 
equation. We have shown that the harmonic oscillator on i/|, as in the other spaces with 
constant curvature, is exactly solvable and belongs to the class of maximally superintegrable 
system. We have proved that all the bounded classical trajectories are closed and periodic. 

The orbits of motion are ellipses or circles for bounded motion and ultraellipses or equidistant 
curve for infinite ones. 

Key words: superintegrable systems; harmonic oscillator; hyperbolic space; Hamilton-Jacobi 
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1 Introduction 

The harmonic oscillator as a distinguished dynamical system plays the fundamental role in 
theoretical and mathematical physics due to many special properties outgoing from its hidden 
symmetry. Together with the Kepler-Coulomb problem they are only one among the central 
potentials for which all classical trajectories are closed (Bertrand theorem) and in quantum me¬ 
chanics all energy state are multiply degenerate (accidental degeneracy). The other consequence 
of hidden symmetry is the existence of additional functionally (in quantum mechanics linearly) 
independent integrals of motion and the phenomena of multiseparability, that is separability of 
variables in Hamilton-Jacobi or Schrodinger equation in more than one orthogonal systems of 
coordinate. It has long been known [9, 12, 26] that the harmonic oscillator problem possesses 
five functionally independent integrals of motion, which generate the separation of variable into 
eight systems of coordinates [11, 17]. In most of them harmonic oscillator admits the exact so¬ 
lution, the fact which makes it attractive to use as a model of molecular, atomic and nuclear 
physics and other branches of theoretical physics. 

The generalization of Kepler-Coulomb system and oscillator problem on the spaces of con¬ 
stant curvature start from the work of Lobachevsky, who first identified the Kepler potential in 
hyperbolic space (two-sheeted hyperboloid) and found the trajectories of classical motion [38] 
(see also the articles [7, 10, 34, 47]). The extension of the harmonic oscillator problem on the 

*This paper is a contribution to the Special Issue on Analytical Mechanics and Differential Geometry in honour 
of Sergio Benenti. The full collection is available at http://www.emis.de/journals/SIGMA/Benenti.html 
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spherical and hyperbolic geometries has already been done in the book of Liebmann [37], who 
also discussed the geometric character of the conics in noneuclidean geometry. The investigation 
of Kepler-Coulomb problem in quantum mechanics was motivated to compare the properties of 
the Coulomb potential in the “open hyperbolic” or “closed” universe to that of an “open but 
flat” universe. Schrodinger [46] was the first who discussed this problem and discovered that 
for “hydrogen atom” on three-dimensional sphere only discrete spectrum exists. Virtually at 
the same time, Infeld and Shild [25] found that in an open hyperbolic universe there is only 
a finite (but very large) number of bound states. The motion in Coulomb field on imaginary 
Lobachevsky space (one-sheet hyperboloid), as shown by Crosche [16], has some peculiarities. 
It is not singularfor any value of variable and their discrete spectra infinite degenerate. The 
essential advance in the theory of systems with hiddensymmetry in the spaces with constant 
curvature was made by Higgs [24], Leemon [36] and Belorussian authors in [2]. They have 
shown that the complete degeneracy of the spectrum of the Coulomb and oscillator problems on 
the three- dimensional sphere and hyperboloid is caused by an additional integrals of motion: 
“curved” Runge-Lenz’s vector (for the Coulomb potential) and Demkov-Fradkin tensor (for the 
oscillator). However, in contrast to the flat space, commutation relations between the compo¬ 
nents of Runge-Lenz’s vectorand Demkov-Fradkin tensor on the sphere and hyperboloid form 
the quadratic or cubic algebra. Later it was proven that these properties are inherent in all class 
of maximally second-order superintegrable systems, which also belong to the Kepler-Coulomb 
and oscillator potentials(see for instance recent review [40] and references therein). 

We recall that in general, in an V-dimensional space, maximal superintegrability means 
that the classical Hamiltonian allows {2N — 1) functionally independent integrals of motion 
(including the Hamiltonian) that are well defined functions on phase space. The first searchof 
superintegrable systems in two- and three-dimensional flat space was done in the pioneering 
works of Winternitz and Smorodinsky with co-authors in [39, 50], later the notion of super¬ 
integrability in the spaces of constant curvature has been introduced in theseries of papers 
[17, 18, 19, 20]. The complete classification of superintegrable systems on the two-dimensional 
complex sphere,which include to real spaces, sphere and hyperboloid, as particular cases have 
done in the work [28]. Some of the superintegrable systems have been constructed on S'^r and Hj\f 
spaces in [23] . We can also mention some articles devoted to the investigationof various aspects 
of both classical and quantum superintegrable systems in the spaces of constant curvature, for 
instance [2, 21, 27, 30, 31, 32]. 

The classical and quantum mechanical systems on the spaces of constant curvature (positive 
and negative) have always drawn a great attention due to their connection with the relativistic 
physics and gravity. The 2D and 3D one-sheeted and SO(2,2) hyperboloids are the models of 
the relativistic spacetime with a constant curvature, namely de Sitter and anti de Sitter spaces, 
which is a crucial point for its wide application in the field theories [42, 49], quantum gravity and 
cosmology [1, 14, 48], integrable Yang-Mills-Higgs equation (or Bogomolny equation) [33, 51]. 
Among other applications we can mention also quantum Hall effect [3] and coherent statequan- 
tization [13]. 

However, as far as we know, the superintegrable systems on imaginary Lobachevski space 
SO(3, l)/SO(2,1), (de Sitter space time dS^^^) on hyperboloid = SO(2, 2)/SO(2,1), 
(Anti de Sitter space time AdS^"''^), have not been studied with the same degree of detail 
and need to be further investigated.lt appears that the first work in this direction (if we do 
not take into account the paper [16]) was the article [8] (see also more general case in [4]) 
where the authors, using the reduction procedure to the free Hamiltonian on the homogeneous 
space SU(2,2)/U(2,1), obtain the eleven different types of maximally superintegrable systems 
on the hyperboloid H^- Later, in paper [22], the superintegrable generalization of harmonic 
oscillator and Kepler-Coulomb potentials covering the six three-dimensional spaces of constant 
curvature (including de Sitter and anti de Sitter spaces) in unified way, parametrized by two 
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contraction parameters defining themetric in each space, have been constructed. In these papers 
the classical superintegrable systems are only identified but have not been solved. Recently, also 
the main properties of two-dimensional harmonic oscillator problem have been investigated in [6] , 
using again two parameters approach, in nine standard two-dimensional Cayley-Klein spaces, 
including the de Sitter and anti de Sitter spaces. 

The present work in a sense can be considered as a continuation of our previous articles [43, 
44, 45], devoted to the investigation of classical and quantum Kepler-Coulomb problem and 
quantum harmonic oscillator problem on the configuration hyperbolic space with constant cur¬ 
vature i?|. The given paper aims to investigate the harmonic oscillator problem on the whole hy¬ 
perbolic space H 2 from the point of view of classical mechanics, which, to our knowledge, has not 
been elucidated in literature so far. This task seems more complicated but also more interesting 
than the analogous problem in the other three-dimensional hyperbolic spaces. It mainly derive 
from the complexity of the space which includes such subspaces as the one- and two-sheeted 
hyperboloids. This study will hopefully also help us to better understand the quantum case. 


2 The hyperbolic space and constants of motion 

A three-dimensional hyperboloid H 2 CJl 2,2 is described by the equation 

4 + 


( 2 . 1 ) 


To be more specific we parametrize the hyperboloid (2.1) using the geodesic pseudo-spherical 
coordinate (r, r, (/?) [29, 43], namely 


zo = TRcoshr, 

Z 2 = R sinh r cosh r cos (p, 


zi = R sinh r sinh r, 

Z 3 = R sinh r cosh r sin tp, 


( 2 . 2 ) 


where r > 0 is the “geodesic radial angle”, r G (— 00 , 00 ), and p G [0, 27r). The connection 
between two sets of coordinates zq ^ —zq corresponds to the complex transformation of radial 
angle r ^ in — r. The system of coordinate (2.2) is valid only for \zo\ > R and the missing 
part of the surface for |zo| < R may also be taken into account if we use another form of the 
pseudo-spherical coordinate 


zq = ±Rcos X, 

Z 2 = R sin y sinh p cos p, 


zi = i? sin X cosh/i, 

Z 3 = R sin X sinh p sin p, 


(2.3) 


where now x £ (“f) f )> ^ (— 00 , 00 ) and p G [0, 27r). It is also easy to see that the two 

pseudo-spherical system of coordinate (2.2) and (2.3) are connected by 




T —)> /i — in/2. 


(2.4) 


Here we shall make use of the pseudo-spherical system of coordinate in form (2.2). To investigate 
the motion in the region \zo\ < R, everywhere below, we will use the transformation (2.4). 

The restriction of the pseudo-euclidean metric = G^^dz^^dz’^G^i, = diag(—1, —1,1,1), 
{p, v = 0,1, 2, 3) on R 2,2 to leads to the following formula 

—^ — sinh^ rdr^ -|- sinh^ r cosh^ ndp'^. 

Then the kinetic energy is given by 

T = — (r^ — sinh^ r (f^ — cosh^ 
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and the canonical momenta can be obtained in a usual way 

dT o dT 2 2 dT 2 2 2 

Pr = = R r, Pt = = —R sinh rf, p^ = = R sinh r cosh rip. 

or OT op 

Thus the free Hamiltonian in the pseudo-spherical phase space {r,T,ip;pr,PT,Pip) with respect 
to the canonical Lie-Poisson brackets 


{f,9} = 


i=l 


_ dg df 

dqi dpi dqi dpi J ’ 


(2.5) 


has the form 


R-iree 


1 


2i?2 


Pr 


sinh^ r 


Pr 


P% 


cosh^ T 


( 2 . 6 ) 


It is clear that isometry group of hyperboloid is given by SO(2, 2) group. The correspon¬ 
ding Lie algebra is six dimensional. The generators of so(2,2) algebra can be written in terms 
of the ambient space R 2,2 coordinates and momenta p^ as 

Cl = -{Z2P3 - Z3P2), C2 = -{ziP3 -h Z3P1), £3 = {zip2 + Z2P1), 

M = {zopi - Zipo), M 2 = -{zqP2 + Z2Pq), Ms = -{zqPs + zspo), (2.7) 


and the Lie-Poisson brackets (2.5) with the help of three-dimensional metric = diag{l,—1,—1} 

reads 


— 9im9jn^mnkCki — 9im9jn^mnkCki {-^ii Cj^ — 9im9jn^mnkMki 

where A: = 1,2,3. There are two Casimir invariants, the first of which vanishes in realiza¬ 
tion (2.7): 

Ci=L-N = N- L = QikMiCk = MiCi — M 2 C 2 — M 3 C 3 = 0, (2.8) 


and the second one is 

C2 = N^ + L‘^, (2.9) 

where 


N^ = N-N = PikMiMk = M^ - Mi - Ml 

£2 = L • L = g^kCiCk = Cl-Cl- £§. (2.10) 

The next step is computing the relationship between the ambient momenta and the geodesic 
polar one. Taking into account that four-dimensional canonical momentum p^ {p = 0,1, 2,3) 

dC 1 


where £ is a kinetic energy in the ambient space R 2 , 2 , we obtain that 
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R.p, = R.^ 


cosh r cosh r cos ip Pr 
cosh r cosh r sin ip pr 


sinh r cos ip sin ip 

sinhr sinhr coshr^"^’ 

sinh r sin ip cos ip 

sinh r sinh r cosh 


Then the generators (2.7) in geodesic pseudo-spherical coordinates and momenta are given by 
the formulas 


A/i = — sinh Tpr + cosh r coth r p ^, 

A/2 = 

A /3 = 

C3 = 


coth r sm ip 

cosh r cos ippr + coth r sinh r cos ippr -\ -^- Pu >, 

cosh r 

coth r cos ip 

cosh r sm ippr + coth r sinh r sm ippr _ _^ _ — Pip , 


coshr 


sm ip 

COSipPr + —^ Pifi, 


coth r' 


cos ip 

L2 =-svaippr - Tp—p^j Ll=p^. 

coth r 


( 2 . 11 ) 


Using now equations (2.9), (2.10) and (2.11) it is easy to see the second Casimir operator C 2 
is related with the free Hamiltonian (2.6) by Ci = —2i?^?7free- Thus all the quantities (2.11) 
Poisson commute with free Hamiltonian (2.6) and are constants of the motion. From the seven 
integrals of the motion {"Hfree,A/",, Tj} only five are functionally independent, because of the 
relation (2.9) and constraint (2.8). Hence the geodesic motion with the Hamiltonian (2.6) turns 
out to be a maximally superintegrable system. 

Let us now consider the spherically symmetric model, namely the Hamiltonian H = 77free + 
V(r), where Hhee is given by equation (2.6) and V(r) is a potential function. It is obvious that 
the Hamilton-Jacobi equation % = £ for any central potential admit separation of variables in 
the pseudo-spherical system of coordinates (2.2) (and (2.3))^. The pseudo-spherical system of 
coordinates corresponds to the subgroup chains SO(2, 2) D SO(2,1) D SO(2). Thus, the central 
symmetry of Hamiltonian % implies the conservation low of the vector L = {C\,C 2 , C 3 ) with the 
scalar product (2.8), which we can interpreted as Lorenzian “angular momentum”. In particular 
the first component of angular momentum Ci = p^p and Casimir invariant of algebra so(2,1): 







together with the Hamiltonian %: 


n 


1 f 2 

2E? sinh^ r 


+ V(r), 


( 2 . 12 ) 


form the mutually Poisson-involutive system of constants of motion. As it follows from the 
equation (2.12): p^/cosh^r — > 0, the quantity in contrast to the motion in Euclidean 

space (or spheres and two-sheeted hyperboloids), can take not only the positive or zero but also 
the negative value. Another difference is that at the fixed values of p^> Lp‘. The existence 
of an additional independent constant of motion C 2 {C 3 then not independent) means that the 
problem is at least once degenerate and the trajectories placed on the two-dimensional surface. 
For the case of positive putting r = 0, or we obtain that the motion takes place 

on the two-dimensional subspace, namely two-sheeted hyperboloid Zq — Z 2 — = B?, while 

for negative L^, we may put p = 0 or = 0 , and restricted to the one-sheeted hyperboloid 
Zq + Zi — Z2 = R^. 

^Beside of the pseudo-spherical system of coordinates (2.2) the Hamilton-Jacobi equation "Hnee = £ and free 
Schrodinger equation on hyperboloid allow the separation of variables additionally in 70th orthogonal systems 
of coordinates (see for details [29]). 
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In the case of \zq\ < R the formulas for so(2, 2) generators (2.11) are changed accordingly to 
the transformation (2.4). We have 




sinm ^ 


Hence by virtue of above relation, the takes only negative value. Without the loss of generality 
we can put (/? = 0 or = 0 and the motion on again restricted to the one-sheeted hyperboloid 
Zq + zf — Z2 = R'^ ■ 

3 Harmonic oscillator potential 

Let us now concentrate on the spherically symmetric model, namely harmonic oscillator system. 
In the article [45] we have extended the Euclidean isotropic harmonic oscillator potential with 
the frequency uj to our space which is given by 


-jy-OSC _ 


fzi + zi- zf 


= 




2 


tanh^ r, jzoj > R, 


tan^ y, \zo\ < R. 


Respectively the Hamiltonian may be expressed as follow 


^OSC _ 


1 


2 R 2 

for l^ol > R, and 


Pr + 


L 2 


sinh 2 r 


, ,2 p2 

UJ XL 2 

H-tann r 


(3.1) 


^^OSC _ _ 


1 


2R? 


(Px + 


L 2 

sin 2 x 


a; 2 R 2 


tan2 y 


(3.2) 


for l^ol < R- 

The Hamiltonian of the harmonic oscillator system, besides the angular momentum L has ad¬ 
ditional integrals of motion quadratic in the momenta, which are associated with the generators 
(A/i,A/ 2 , A/ 3 ), the so called Demkov-Fradkin tensor [9, 12 ]: 


K Zq 


T^ik = T^kh i,fe = l,2,3. 


The components of Vik tensor Poisson commute with Hamiltonian of harmonic oscillator (3.1) 
and (3.2), but not necessarily with each other. In the pseudo-spherical coordinates the diagonal 
components of this tensor has the form 

^11 = sinh2 r tanh2 r, V 22 = + oJ^R? cosh2 r cos2 tp tanh2 r, 

R^ R^ 

^33 = cosh2 r sin2 ip tanh2 r, 

R^ 

so the harmonic oscillator Hamiltonian is given by 

L 2 


2 R 2 ' 


>^080 


Vii -I- V22 + iRss 


(3.3) 

















Harmonic Oscillator on the SO(2,2) Hyperboloid 


7 


In addition to this, the Demkov-Fradkin tensor has the algebraic properties 

Y^CiVik = Y,'^k^C-i = ^, /s = l,2,3. (3.4) 

i i 

It is clear that the ten integrals of motion {%,Ci,T)ik} cannot be functionally independent 
because of the relations (3.3) and (3.4), and that 


{CiDii] = {£ 27 ^ 22 } = = 0 . 


Only five integrals of motion, which we can choose as {77, £ 1 , £ 2 , are functionally in¬ 

dependent. Thus 77°®'^ is a maximally superintegrable Hamiltonian. The components of angular 
momentum and Demkov-Fradkin tensor forms the quadratic algebra. The nonvanishing Poisson 
brackets have been presented in Appendix A. 

In the contraction limit 7?—>-oo the hyperbolic space turns into the Minkowski space 

Let us pass to Beltrami coordinates 


Xi = R— = R 
zq 


Zi 


y/R? z\ 


Zrzl-z\ 


i = 1,2,3. 


(3.5) 


Then, at the limit 72 —t- 00 we have that 


lim H“Dr) 
R^oo 


^ ( 2 , 2 I 

y(-Xi +X2 +X3), 


which can be interpreted as a harmonic oscillator potential on the Minkowski space 

(X1,X2,X3). 


4 Integration of the Hamilton—Jacobi equation 


The Hamilton-Jacobi equation, associated with the Hamiltonian (3.1), is obtained after the 
substitution —)• dSjdni, where /r* = {r,T,ip). Therefore we get 


77 = 


2722 


dr 


1 


+ 


1 


sinh^ r \ dr J sinh^ r cosh2 r V 


dS 


2 7?2 


> + 


uj^R 


tanh^ r = E. 


This equation is completely separable, and the coordinate ip is cyclic. We look the solution for 
the classical action S{r,T,p) in form 

5(r, r, (p) = p^p J- 5i(r) + S 2 {t) - Et, 

and obtain 


dS2 

dr 


P% 


- 


cosh^ r 
2 


= -L^ 


1 fdSiY uj^R'^ ,2 

-|-tanm r -|- 


£2 


2722 \ dr J ' 2 ' 2722 sinh^: 


= E. 


The “quasi-radial” equation (4.2) describes the motion in field of effective potential 


2 d 2 


t/eff(r) = 


oj^R 


■ tanh^ r + 


£2 


2 2722 sinh^' 


(4.1) 

(4.2) 

(4.3) 


At the large r ~ 00 the effective potential ?7eff(r) tends to a constant value equal to a;^72^/2, 
whereas the behavior at the point r = 0 is determined by the angular momentum £2. 
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Figure 1. Effective potential C/eff(r) in case of 0 < for value of = 0,1/16,1/8,1/4; 

uj = R=l. 


l’eS(r) 



0 1 2 3 4 5 

Figure 2. Effective potential Uef{{r) in case of > uP'R'^ for value of = 2, 3,4; a; = i? = 1. 


In case 0 < potential (4.3) has a minimum at tq = tanh ^ ^ 1? juP'Ti^ (see Fig. 1), 

and at this point 

0 < C/eff(ro) = cuVP - ^ (4.4) 

where equality is possible only in case of I? = 0. For > uP'R^ the potential ?7eff(?’) is 
repulsive on the whole semi-axis r G [0, oo) (see Fig. 2). In the case of negative 1? the effective 
potential (4.4) is attractive and has a singularity for a small r as ~ (see Fig. 3). 

For the region l^ol < R the differential equations (4.1) and (4.2) are transformed to the 
following ones 


v% 

\ dfj, J sinh^ fj. 

1 a;2R2 


-L^ 


tan^ X -|- 


2i?2 gin^ X 


-E. 


The first equation admits only negative value of L^. Therefore we take into account the motion 
inside the region \zo\ < R when investigate we the case of negative value of L^. 
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Figure 3. Effective potential Uef{{r) in case of < 0 for value of = —1, —2, —3; uj = R = 1. 


Integrating now equations (4.1) and (4.2) we get 


Si{r) 

52 (r) 



- tanh^ 


pI 


cosh^ r 


dr. 


r 


L 2 

sinh^ r 


dr, 


(4.5) 

(4.6) 


Since we are interested only the trajectories we will follow the usual procedures [35] and consider 
the equations 


dS^_dSi_ __ dS__dS^ ^-r 

dE dE ^ dE^ dE^ dE^ ’ dp,p ^ dp^ 


(4.7) 


where Iq, ipQ and (5 are the constants. 


4.1 Integration of quasi-radial part 

From equations (4.5) and (4.7) we get that 


1 

t - to = - 

OJ 


tanhrdr 

sj- tanh^ r + 2 {EluEB? + E^ tanh^ r - E^ 


(4.8) 


Below we consider separately all four cases: 0 < E^ < uER^, > uER^, L? < 0 and = 0. 

1. The case 0 < < u'^R^. For the roots in the radical expression of denominator in (4.8) 

we have 


Xi,2 


{2R^E + E^) ± ^{2R^E + L2)2 _ AL'^u'^R^ 
2^2 


(4.9) 


where X = tanh^ r G [0,1]. It’s obvious that the radicand in equation (4.9) is positive for 
any values of energy E > F'min = t4ff(?’o) and equal zero for E = F^min- Thus the roots Xi _2 
(Xi < X 2 } are positive. It is easy to see that for F^min < E < u'^R? j2 both roots satisfy the 
inequality condition 0 < Xi < X 2 < 1. At F^ > uElEI2\ 0 < Xi < 1 < X 2 and equality ^2 = 1 
is possible only for E = oj'^R?I2. The bounded motion exists exclusively for F^min < E < uER?I2. 
Below we will consider separately all possible cases, namely: Smin < E < oj'^R?/2, E = Smin, 
E > u?R^l2 and E = u?R^l2. 
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A. Performing the integration in formula (4.8) we get for Pimin < E < j2 
2 uj‘^R^ sinh^ r = (1 - 2 E/uj'^R^)-'^^{2E - L^/R^) 

+ yJ~{2E^^E^jR^f^^AL^ sin [2ujy^l - 2EliJ^R?[t - to)] }• 
Thus the motion is bounded and periodic. The period is given by 


vr 


T{R) = - 


w ^1 - 2Elu:‘^R? ■ 


(4.10) 


The total frequency loq = wy^l — 2EjuP'R? and unlike the motion in Euclidean space, depends 
on the energy of particle E and curvature of the space k = —1/i?^ as a parameter, but it is 
constant for each of the orbits at a fixed value of the energy^. This property is common to all 
closed orbits of superintegrable systems on the spaces with constant curvature. The contraction 
limit i? —)■ oo give us the correct Euclidean period: T{R)ji^oo = 5- period of motion on 
always larger than in Euclidean space by the factor: 1/y^l — 2E/uj‘^R? and tends to infinity at 
the limit E —)• uj'^R'^ j2, that is the closed orbits changes to the infinite open ones. 

B. In the case of minimum energy: E = Emin = t^eff(^’o) or E min = ujVi? — L^/2E^ the 
integral in (4.8) is not defined and we must solve directly the equation (4.2). Erom equation (4.2) 
we obtain 


(dSi 
V dr 


2 


E^ coth r — Vtanh > 0, 


or dSi/dr = 0 and tanh^r = ^/iXJuER^. Therefore 


r = tanh 


-1 


1- A 1- 


2E 

I ’ 


(4.11) 


i.e., the trajectories are circles. Here from two values of V~E allowed by equation E = ?7eff(?'o)) 
we choose the smaller one VlE = ojR? ^1 — y^l — 2E/uj'^R?^ because it satisfies the condition 

0 < < uj'^R^. In case of contraction limit E —)■ oo we obtain E = E min = ujVlZ and 

r = y/E/oj. 

C. In case of E > oj'^R?/2 after integration in (4.8) we have 
2oj‘^R^sm\Er = {2 EIu?EZ - l)~^\^[L^/R^ - 2E) 

+ ^( 2 E + E2/E2)^ -4L2a;2 cosh [2uj^j2E/uj‘^R^ - l(to - i)] }, (4.12) 

i.e., the motion is not bounded. 

D. For the limiting case of E = uER?I2 the roots of denominator are Xi = L? juERZ^ X 2 = 1, 
thus /uErZ < tanh^ r < 1 and motion is not bounded because of tanh“^(L^/a;^E'^) < r < 00 . 
The integration in (4.8) yield 

cosh^ r = (1 — ^ (l — E^/cj^E^)(t — to)^- (4-13) 

2. Let us consider now the case of E^ > cj^E^ (see Fig. 2). From equation (4.8) we get 
that the only possible value for energy is E > a;^E^/2 and the roots satisfy the inequality 

^The Euclidean harmonic oscillator is a classical example of an isochronous system [5]. The period of motion 
of Euclidean oscillator depends only from frequency and is the same for all orbit. 
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0 < Xi < 1 < X 2 . Thereby, the equation of motion is determined by the formula (4.12). The 
motion of particle is limited only by the point r m\n = tanh“^ V^T, i-e., it has the ability to go 
to infinity. 

3. Let us consider finally the case of < 0. From the equation (4.8) we have that the roots 
of denominator are 


^ 1,2 = 


{2ER‘^ - |L2|) ± {2ER^ - \L‘^\f + 4\L‘^\u;^R^ 


2^2 


where again X = tanh^ r G [0,1]. It can be seen that Xi < 0 < X 2 is independent of the value 
of A and energy E. For the region E > lo'^R'^ j2 one of the roots is X 2 > 1, so the radicand is 
positive for any values of variable r, including the point r = 0: r G [ 0 , 00 ). The same situation 
develops for region E < oj‘^R‘^/2, where r G [0, tanh“^\/X 2 ]. Therefore in case of negative A the 
particle can penetrate from the region zq > R io < zq < R. 

Performing the integration in formula (4.8), we have for E < uP'R? j2 


sinh^ r = 


+ 


while for E > uj^R^ j2 


2 REE^\L^\ 

2 i?2(t^2^2 _ 2E) 

y/{2R^E - |L2|)2 + 4|L2|a;2i24 


2 i?2(t^2^2 _ 2E) 


sm 


2u;y^1 — 2E/uj‘^R?{t — to) 


(4.14) 


sinh^ r = 


2R^E + \L^\ 

2 i?2(t^2^2 _ 2E) 

y^{2R^E - |L2|)2 + 4|L2|a;2i24 


+ 


2R‘^{2E - u}‘^R‘^) 


cosh 


2 u:'j 2 E juP'R? — l(t — to) 


(4.15) 


From the formula (4.14) it follows that the motion E < (jj‘^R? j2 is bounded and periodic with 
period (4.10). Below we will construct the bounded trajectories lying on the whole hyperboloid, 
namely not only in the region \zo\ > R, but also \zo\ < R- In case when E = oj'^R? the integration 
in (4.8) leads, up to a transformation —)• —\L?‘\, to a result similar to the formula (4.13). 

In the limiting case of = 0 the formulas (4.14), (4.15) and (4.13) are simplified. For 
0 < i? < zj'^R?/2 we get 


sinh^ r = 


2 e >2 


1 - (“Vl - - to) - , 

while in case of E > tJ^R? j2 

sinhr = sinh {uj2 E/R^ - l(to - t)) • 

Finally for E = u'^R?/2 we obtain sinhr = uj{t — to)- 

4.2 Integration of the angular parts 

1. Let us first consider the case when > 0. From (4.5) and (4.6) we obtain 


dSi 


dr 


sinh^ r-y^ 2R?‘E — uj'^R'^ tanh^ r — j sinh^ r 


(4.16) 
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dSo 


dr 


2 J ^_L2+p2/cosh2r 

The integrals can be easily calculated to give [15] 

dS2 


(4.17) 


1 

dSi 1 


arcsin 


sinhr 


5L2 4^ 

Here we require 


arcsin 


2 L 2 coth^ r - {2ER^ + L^) 
y/{2ER^ + L2)2 _ 4L^co‘^Ri 


-\ pI/L'^ - 1 < sinhr < Jp%/L‘^ - 1 


and 


\2L^ coth^ r - {2ER^ + L^) | < W [2ER'^ + L"^) - AL^lo^R-^. 


(4.18) 


The condition (4.18) is equivalent to zi < cothr < Z 2 , where 2 : 1,2 are the roots of denominator 
in integral (4.16): 


Z 12 = 


{2ER^ + L2) ± ^{2ER‘^ + L2) - AL^uo^R^ 


E > E^in = ioVl^ - L^/2R^. 


The final condition Z 2 > I implies that and E > ui‘^R^j2 or 0 < < uj'^R'^ and 

E i> P/min- 

Therefore for dS/dL^ we have 


1 


dS 


aL2 4^/])2 



2 L 2 coth^ r - {2ER‘^ + L^) 


{2ER^ + L2)2 _ 4L2a;2^4 


— 2 arcsin 


sinhr 


pI/l^ -1 


= /?■ 


(4.19) 


Next, from (4.6) and (4.7) we obtain 


dS _ p^dr 

^ cosh^ ry'^— 

and hence 

tanhr = ^1 - E^/p^^s\n{tpQ - ip). 


tanhr 

= = </? + arcsin ■— — = = ipQ, 




(4.20) 


(4.21) 


2. Let us consider the integration in formulas (4.16), (4.17) and (4.20) in the case < 0. 
Instead of equation (4.19) we obtain [15] 



2lL2|coth2r + {2ER^ - Jy!]) 

5 L 2 4 ^ arcco.il 

_^J{2ER^ - 1 L 2|)2 ^ 4|^2|^2^4^ 
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— 2 arcsinh 


sinhr 




= / 3 , 


and 


sin((^o — ‘f) = — _P^£__ 


(4.22) 


(4.23) 


with the restriction for r: 


coth r > I - 


ER^\ [7l ER‘^y 


The limiting case of = 0 could be easily calculated directly from equations (4.22) and (4.23). 
So, we get 


dS \/ 2E coth^ r — uj'^R? sinh r 

< 9^2 ^ 2=0 

with the obvious restriction coth^r > oj'^R^I2E. 


sinhr = tan((/7o — 


(4.24) 


5 The trajectories for > 0 

From (4.19) and (4.21) we have 


,2 I'ER^ 1\ UER'^ iV . /T9n\ 

coth^ r = 1 + - I + \ l i + - ) - fY~ + 4:^1?j3 ), 


L 2 2 


L 2 


where 


(5.1) 


= arcsin 


sinhr 


- 1 ^ 


= arcsin 


^1 + L‘^/pI cot2((^o -<p)) 


Now we can rewrite the equation (5.1) in form 
tanh^ r=— ^ 




I 1 
L 


(5.2) 


(5.3) 


Thus we see from (5.2) that the dependence of angle r in the equation of trajectories (5.3) can 
be eliminated. On the other hand from the formula (4.21) it follows that the motion of particle 
on the hyperboloid is restricted to the additional condition 


zi tanh r 
Z 3 sin (p 


y^l -L2/p2. 


Therefore, without the loss of generality we can choose t = 0 ov L'^ = p^. Taking into account 
that the formula (5.3) is invariant about transformation r ^ in — r we can conclude that all 
trajectories of motion, given by this formula, lie on the upper (zq > R) or lower {zq < —R) 
sheets of the two-sheeted hyperboloid: Zq — Z 2 — = R?. Obviously they are symmetric with 

respect to transformation zq —>■ —zq. 
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Putting now in (4.21) we obtain that V' = ifo ~ V’) the formula (5.3) gain the 

following form (equation of orbits) 


tanm r = 


P 


1 + e{R) cos 2 (/ 3 ’ 
where we use the notations 


(5.4) 


piR) = 


fER^ 1 
2 


> 0 , 


e{R) = Jl- 


{2ER^ + L‘^Y 


< 1 , 


(5.5) 


and choose (pQ = —2^\fA + ^ that the points (/? = 0 will be the nearest to the center. It is clear 
that radicand is always positive because of PI > ?7eff(?’o) for 0 < A < uP'R^ and E > u'^R?/2 for 
A > cj^R^. 

It is well-known that as in the Euclidean plane it is possible to introduce the conic (section) 
on the two-dimensional spaces of constant curvature [7, 10, 34] (see also the definition of curves 
on the two dimensional hyperboloid in [41]). The conics on the spaces with constant curvatures 
are the curves of the intersection between two-sheeted hyperboloid (or sphere) and second order 
quadric cone with the origin in the center of hyperboloid (sphere). Geometrically the conic 
on the spaces of constant curvature possesses many properties characteristic of conic section in 
Euclidean plane, particularly we can speak about the focuses E^ and E 2 and can determine the 
conic as the point set, from which the sum (ellipses) or difference (hyperbolas) 2 a of distances ri 
and r 2 to two given points (focuses Ei and E 2 ) are constant. 

Let us now analysis of the oscillator orbit (5.4). The formula of trajectories 
written in more convenient form 


(5. 

or in term of the Beltrami coordinate (3.5): 


1 cos^ cp sin^ p 

tanh^ r ' 


(5.4) may be 


Xn 


_^ I _^ _ DZ 


(5.7) 


where the constant A and B are 


„2 ^ P{R) a2 ^ P(.R) 

l + e(R)’ l-e(R)’ 


0 < < A^. 


(5.8) 


The orbit equation of the type (5.6) has been studied in detail in the paper [ 6 ] (see also [10]) 
at the investigation of two-dimensional harmonic oscillator in the space of constant curvature 
in polar coordinates. The curves (5.6) are always conic on the hyperbolic plane, but its type 
depends on the value of A and B. It is obvious that if the value A^ > 1 and B^ > 1, then 
for any polar angle p it follows that tanhr > 1, and this case cannot produce any oscillator 
orbit. In the case of B^ < A^ < 1 the conic (5.6) takes the form of hyperbolic ellipses. The 
quantities A and B are related to the lengths of the large and small semiaxes a and b, running 
the interval [ 0 , 00 ), defined as the values of r at p = 7 r /2 and p = 0. Then the values A, B 
can be written in term of hyperbolic tangent of a, b: A"^ = tanh^ a and B^ = tanh^ b and the 
equation of orbit (5.6) is 


1 

tanh^ r 


cos^ p 
tanh^ b 


sin^ p 
tanh^ a 


(5.9) 
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In the contraction limit i? —?• oo we have r ^ f/R where f = is the radial variable in 

the Euclidean plane. Taking into account the limit 


eiR) 


e = \ 1- 




R?p{R) —)• p = —, 


E2 ’ ” ' "" E 

we get from (5.7) that the equation of trajectories transforms into the oscillator one 


Xn 


= 1 , 

H2 H2 


B^= ^ 


1 + e’ 




P 


1 — e 


The next interesting case is when < 1 < A^. This conic (5.6) is neither the ellipse nor the 
hyperbola. Following the paper [ 6 ] we will call this conic as the ultraellipse. Only one semi¬ 
axis b belongs to the hyperbolic plane and the next one formally is not on the real distance. It is 
possible to introduce a new “semiaxis” a (situated on the complex plane on the line a = a+iir/ 2 ) 
which related with the quantity Ahy = cotha. Thus, instead of (5.9) we have the conic 


I 


COS^ P 1 2 2 

-^ = - 2 - 1 " tanh asm ip. 

tanh r tanh b 


(5.10) 


There is a joint point of two conics (5.9) and (5.10), namely ^.2 = I (a —?• oo or a —?• oo). In this 
case the conic is given by 


I 


COS^ P . 9 
^ + sm^ ip. 


tanh^ r tanh^ b 

This conic is an equidistant curve with equidistance b from the axis Z 2 [ 6 ] . 

Let us now consider all the possible trajectories of motion depending on the energy and 
angular momentum A. 

A. First we consider the case when UeiAo) < E < Ar'^ j2 and 0 < < AR'^. It is clear 

that 


-1 


B^ < A^ = 


'ER^ , 1 , 

IX 2 J “ V [jl? 


fEA^ ly w2i?4 

^ 2 ) ~ ~1X 


< 1 


and the oscillator orbits are described by the equation (5.9). Denote the minimum b = r^\^, 
(p = 0 ) and maximum a = rmax, {p = 7 r/ 2 ) points on the orbit as a distance from the center of 
field. From (5.8) and (5.9) we have 


tanh^ Tmin = 


p 


I + £{R) 
and correspondingly 


tanh^ Tmax = 


P 


I -e(i?)’ 


= COth ^ 


fER^ I 




' ( EB? 1 


+ 2j+VV^ + 2 


W2i?4 

XX 


, 


= coth ^ < 


fER^ I\ UER"^ iV AR^ 




Thus we find that the trajectories of motion are ellipses lying symmetrically to the point zq = R, 
zi = Z 2 = Z 3 = 0 on the upper sheet of the two-sheeted hyperboloid (see Fig. 4). 
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Figure 4. The figure shows the elliptic trajectories lying on the upper sheet of the two-sheeted hyper¬ 
boloid Zq — Z 2 — zl = zq > R for the value £ = 0.3 and p = 0.3,0.4, 0.5. 



Figure 5. The cyclic orbits: £ = 0 and p = 0.2,0.5, 0.8. 


B. In case of minimum energy E = E m\n = fleff(^o) we have from (5.5) that e = 0 and 

p = and consequently tanh^ r = = uA/y/l?. Thus the orbits are circles 

with the radius given by the formula (4.11) (see Fig. 5). 

C. For the case of energy values E = we get that 


p{R) = 


2A 


+ A ' 


e{R) 


AR^ - A\ 
AR^ + A ’ 


therefore for 0 < A < AR'^ we get = A/AR^ < 1 and = 1. The conic is 


1 


2d4 


AR 


tanh^ r A 


cos^ (f -)- sin^ (f, 


which represents the equidistant curves (see Fig. 6). The minimal distance rmin from the center 
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Figure 6 . The figure shows the equidistant orbits lying on the upper sheet of the two-sheeted hyperboloid 
• 2^0 ~ -^2 ~ -^3 = zq > R with the value of pairs (p, e) = (1/3,2/3); (2/3,1/3); (8/9,1/9). 


is given by the formula 


^min 


coth ^ 



Let Then = 1 and the conic is a “largest” circle with radius r = oo. For 

the case we obtain that A? = > 1. Then from the formula (5.6) it 

follows that tanhr > 1 and no any oscillator orbits exist. 

D. For the energy E > j2 it is easy to see that for any positive > 0 






2 ) L 2 I 


> 1 , 


< 1 . 


The motion of a particle is determined by the equation (5.10) where tanh^ d = \/A?. The 
trajectories are ultraellipses and describe the motion of a particle from the minimum point 

ER^ 1 \ I fER^ iV 

to infinity (see Fig. 7). On the other hand side B^ ■ = uAR^ jl?^ so that for I? < uP'R^ 

we get XjA? < < 1 , whereas for > ijAR^ follows that B^ < XjA? < 1 and the value of 

= uAR^ or B^ = \jA^ separates two set of ultraellipses. 

Let us also note that the in contraction limit i? —)■ 00 these orbits corresponds to the Euclidean 
oscillator orbits with the large values of energy (the straight line = B'^). 


= coth 


-1 


6 The trajectories for < 0 

To simplify further formulas we set first = 0. Then, from equation (4.23) it follows that 
the motion occurs at a constant value of the azimutal angle p = ipo that is limited by the 
condition Z 3 /Z 2 = tan</?o. To further simplify it is enough to choose po = 0 01 c = tt . Thus we 
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Figure 7. The figure shows the ultraellipses lying on the upper sheet of the two-sheeted hyperboloid 
Zq — Z 2 — zl = , zo > R for the value e = 0.8 and p = 0.2,0.5, 0.8. 

get that trajectory of the motion lies on the one-sheeted hyperboloid Zq + zf — Z 2 = R?. The 
formula (4.22) gives us the equation of the trajectory in the region zq > R: 

,2 A ER^\ [7l ER^V a;2R4 , 

’■=(2"iiyj + v(2“Tyrj (6,1) 

Performing the further transformation r —)> zy and r —^ — z 7 r /2 in formula (6.1), we obtain the 
equation of the trajectory in the region 0 < zq < R: 

2 /I ER^\ [7l ER'^V N 

cot x = + ^^cosh( 2 /z + 4vlA/3). ( 6 . 2 ) 

In the formula of trajectory (6.1) we must distinguish two cases, namely for the value of energy 
E < oj'^R? jl and E > lj‘^R^ j2. 

In the first case E < uj‘^R‘^ j2 from equation (6.1) it follows that for any value of the variable 
r G (— 00 , 00 ) we have that cothr > 1. Therefore, the trajectory of the motion extends from the 
point r = 0 at the r —)• —00 (zq = R, zi < 0, Z 2 > 0) to its maximum 


^max 



ER^\ [7l P;R 2 \ 2 ^2^4 


at the point r = —2y^\L‘^\f3 and then goes back to the point r = 0 when t ^ 00 {zq = R, zi > 0, 
Z 2 > 0). Further on, the particle penetrates through the point zq = R from the region zq > Rio 
the region 0 < zq < R, which, as it follows from the equation (6.2), corresponds to the value of 
angles /x —>■ 00 and y —>■ 0, (zq < R, zi > 0, Z 2 >0). Further trajectory extends to the maximal 
value Xmax: 


-1 




ER‘^\ 

W\) 




uj^R^ vr 
|L2| - 2’ 


Xmax — cot 


2 |L2| 
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Figure 8 . The trajectories of motion in the case of |L^| = 1] E = —3/2, —1/2,1/4,1/2,3/2; uj = R = \. 

at the point /r = and then continue to — )• —oo, y —)• 0 (zq < i?, > 0, Z 2 < 0). 

After, the particle again passes the point zq = R and penetrates to the region zq > R. Further 
using similar reasoning it can be shown that the trajectories in case of -E < lo'^R‘^/2, are a closed 
curve lying on the one-sheeted hyperboloid Zq+Zi—Z2 = R? , zq > 0 , so the motions are bounded 
and periodic. The same situation takes place for the case of zq < 0. 

In the case of E > u'^R'^ j2 it is easy to see that the inequality 


2 ' “[Af J ^ |L2| - 2 ~\IJ\ 

is valid. Thus the trajectory of the motion, depending on the sign of variable r is split into two 
paths. One of the paths begins from the large r at the minimal point 


1 


'Tmin = -2^/\lJ\l3 + - cosh ^ 


1 I 

2 + Wf 


1 _ ER'^ 

2 -\IE\ 


+ 


W\ 


and continues to the point r = 0 at r —)• oo (zq = .Ej -^2 > 0). Then the trajectory passing the 
part of 0 < goes back from (zq = R, Z 2 < 0) at the point r = 0, r ~ oo to r G oo at Tmin- 

The second path is symmetric with respect to axis zi. Thus the trajectories of motion in the 
case of E > oj'^R?/2 are not bounded. Some examples of trajectories for the fixed negative 
and various values of energy E, are presented on the Fig. 8 . 

In the case L?' = f) if is easy to get from (4.24) 


cotm r = 


a;2i?2 

2E 


+ Ri/E {2/3 — tan (f/pip)^ 


with ifQ = 0. In the case of E < uj'^R'^ /2 the bounded motion takes place r min = 0 {^p = e/2) 
and Tmax = coth~^ fJ ‘^ 2 E^ ~ arctan2/lp^), whereas for the E > u‘^R?/2 the orbits are 

infinite: r G [0, oo). The trajectories of the motion can be presented on the hyperbolic cylinder 
Zq — Z 2 = R?, zf = Zq, Zq > R (see Fig. 9). 
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Figure 9. The bounded and infinite trajectories of the motion for = 0 lying on the hyperbolic cylinder 
^0 ~ ^2 = ^0 > R- The figure shows the cases E = 0.2,0.5,0.8; w = i? = = 1. 


7 Conclusion 


We have shown that the notion of harmonic oscillator problem can be extended not only to 
the sphere and two-sheeted hyperboloid but also to the hyperbolic space It was proved 
that the harmonic oscillator problem on is exactly solvable and also belongs to the class 
of superintegrable systems. We have constructed the dynamical algebra of symmetry for this 
system, which is nonlinear and quadratic (so-called Higgs algebra). We completely solved the 
Hamilton-Jacobi equation for harmonic oscillator problem in the geodesic pseudo-spherical sys¬ 
tems of coordinates. It was shown that for positive value of the Lorentzian momentum > 0 all 
trajectories of motion lie on the upper (or lower) sheets of two dimensional two-sheeted hyper¬ 
boloid Zq — Z 2 — z'^ = R^. These trajectories are always conics centered in the origin of potential 
r = 0. For the special values of energy Fimin < E < uP'R?(2 and momentum L? < uj'^R^ all the 
orbits are ellipses (or circles for E = Fimin)- In case when E > iJ^R? j2 independently of the 
value of I?, the oscillator orbits are ultraellipses or equidistant curves ioi E = uP'R?/2. We have 
seen that in case of negative values of Loreinzian momentum F^ < 0 the oscillator orbits lie on 
the one-sheeted hyperboloid Zq + z1 — Z 2 = R^ and are bounded and periodic for E < uj‘^R? /2 
and infinite for E > The similar situation is valid for = 0, but in this case the 

orbits lie on the hyperbolic cylinder Zq — Z 2 = R'^, zf = z^. 

Let us make short comments concerning the connection of the classical and quantum case. 
The quantum-mechanical counterpart of the angular momentum operator (2.7) comes through 
the replacement —)■ —id/dZfj, and is given by 

Li =-i{z 2 d 3 - Z 3 d 2 ), L2 =-i{zid3 + Z3di), L3 = i{zid2 + Z2di). 


Then in the pseudo-spherical coordinates (2.2) the operator takes the form 

19,5 1 92 

coshr-r- 


P = Ll-p2-Ll = 


cosh r 9 t 


9t cosh^ r dp , 

and coincide with the Casimir operator of SO(2,1) group. Thus the Schrodinger equation for 
the harmonic oscillator potential can be written as 


1 9 ,2 

„ — smh r— + 
smh^ r or Or 


2R^E- 


P 


sinh^ r 


— ipRl^ tanh^ r 


T = 0, 


(7.1) 
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and solved by separation of variables via the ansatz 'if{r,T,ip) = TZ{r)y {t, ip) . The pseudo- 
spherical function T is a eigenfunction of operator L^y = i{l+l)y which describes the quantum 
geodesic motion on the two-dimensional one-sheeted hyperboloid. The spectrum of i can take 
as well as the real values: i = 0,1,... (discrete series of representation of SO(2,1) group) 
and complex value t = —1/2 + ip, p > {) (continuous principal series). In the first case the 
eigenvalue of Lp' operator is positive and in the second one negative. The exact solution of the 
Schrodinger equation (7.1) for the positive eigenvalues of operator tP has been constructed in 
the previous paper [45]. It was shown that as in the caseof two-sheeted hyperboloid, the energy 
spectrum contains the scattering states and a finite number of degenerate bound states. This 
fact coincides with the existence of closed and infinite orbits for positive Lp in classical case. 
We have not considered in the article [45] the quantum motion in the case of negative eigenvalue 
of because of the strong singularity at the center of harmonic oscillator potential, although 
it is clear that the system has a discrete spectrum. This work is in progress. 

Finally, we wish to emphasize that the Kepler-Coulomb and harmonic oscillator potentials 
are the “building block” upon which most of superintegrable potentials can be constructed. 
Thus the investigation of these systems is important for the further study and understanding of 
more complicated superintegrable systems in the hyperbolic space 77|. 


A Symmetry algebra 

The nonvanishing Poisson brackets between the components 
and Cp. 


{T> 12 , Cl] = —T>13, 

{7?13j -Cl} = T’i 2 , 

{7^23, -Cij = 7^22 — 7 I 33 
{Pii, £ 2 } = —27?i3, 
{7?22, -C3] = —27?i2, 

The same between 

{7?ii 
{7?ii 
{7^22 
{7^22 
{C’as 
{C’as 
{V 12 
{V 12 
{7?13 


{7^12, £2} = —£’23) 

{7^13) £2} = —7?ll — 7?33 
{7^23) £ 2 } = -7^12: 
{T’lijCsj = —2 T’i 2, 
{£’33, £1] = 27?23, 


of Demkov-Fradkin tensor Vij 

{ 7 ?i 2 , £ 3 } = -Till — D 22 , 

{7^13) £3}] = -7^23) 

{7123) £3} = -7^13) 

{'D 22 , £1} = —27?23) 

(7^33) £2} = -27?i3, 


7?12} 

£’23} 

7?12} 

7^23} 

7?12} 

7^23} 

V 13 } 

7 ^ 23 } 

7 ^ 23 } 


= 2 uj £3-1- ^£37^1!) 

= ^(£ 2 £’ i 2 + C3V13), 

= 2 uP C 3 — ^£ 37 ^ 22 ) 

= 2lXp Cl — j^Ci'D22, 

= -^(£27123 - £i7?13)) 

= —2a;^£i -|- j^CiV^s, 


= -{2uj^ - 


4R^ 


+ 


= [2uj^ - 




= -\ 2uj^ - 


AR^ 


£3 + 


{£> 11 , £> 13 } — 2uP C 2 + j^C2'Dii, 

^(£ 37^23 + £i7Ii2) 

^£i£’23, 

{£’33)7Ii3} = C 2 — j^C2T>33 
{'^33,'^n} = -^£27^13) 

(£i 7 ?ii -I- £27)12 + £37)13)) 

(£i7)i2 + £27)22 — £37)23) ) 


^ ( —£i 7 )i 3 -I- C 2 V 23 — C 3 V 33 ) 


{7)11,7)22} 
{7)22) 7)13} = 
{7)22)7)33} = 
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